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Abstract 

Under the sub-wavelength condition (w < A/2), an analytical model of light 
transmission and reflection through ID metallo-dielectric gratings of real 
metals has been developed. It has been shown that the transmission intensity 
associated with the Fabry-Perot (FP) resonance of a ID metallo-dielectric 
grating of a real metal decreases with the increasing grating thickness and 
the dielectric constant of the ridge material. Further, it has also been demon¬ 
strated that the intensity of the FP resonance increases with the increasing 
slit width while it is independent of the grating period (P) and the incidence 
angle (when P « A). 
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1. Introduction 


After the discovery of Extraordinary Optical Transmission (EOT) through 
ID metallo-dielectric gratings, a number of efforts, both numerical and semi- 
analytical, has been undertaken to understand the physics behind this excep¬ 


tional phenomenon 


-9]. Numerical approaches of modeling EOT devices 


include the transfer matrix formalism jl, 5, 10, 11], the Finite Difference 
Time Domain (FDTD) method 


sis (RCWA) 


3, |S 




12], the Rigorous Coupled Wave Analy- 


8|, 13], and the Exact Modal Analysis (EMA) |1| . Both 
the transfer matrix and FDTD methods solve Maxwell’s equations directly 
using grids and produce accurate results. Other features of these two meth¬ 
ods include the ability of visualizing electric and magnetic fields in a plane, 
both inside and outside of a structure understudy. Further, flow of energy 
(Poynting vector) at any location at a particular instant of time can also 
be calculated and monitored. However, the transfer matrix formalism and 
the FDTD method are computationally very demanding 15]. In contrast, 
the RCWA and the EMA are the two closely related numerical techniques 
which are widely used to analyze ID structures { 2 , s|, 5, 6, 8,16]. The RCWA 
expresses the dielectric function of the overall grating structure and the Elec¬ 
tromagnetic (EM) held inside the grating region (—h < z < 0, see Fig. [T|) 
in terms of Fourier series jj0. With the help of Maxwell’s equations a space 
harmonic of the EM held can be expressed as a function of the remaining 
harmonics of the EM held and harmonics of the dielectric function. After ap¬ 
plying proper boundary conditions at z = 0 and z = —h interfaces, rehection 
and transmission efficiencies of various diffraction orders can be calculated 
[l3 ]. I 11 the case of EMA, electromagnetic helds inside the grating region are 
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expressed as waveguide modes. Profiles and wave numbers of various modes 
are determined numerically by solving a transcendental equation obtained 
using the stratified medium theory 14]. Once profiles and wave numbers 
of different modes are available, Maxwell’s equations in combination with 
boundary conditions (those like the RCWA) are solved to find reflection and 
transmission coefficients. At this point it is relevant to mention that al¬ 
though the numerical techniques discussed here produce good results, they 
are devoid of physical insights. Specifically, relationships among transmission 
and reflection coefficients, and w, P and h are not clearly understood. Fur¬ 
ther, simulation techniques discussed above require significant computational 
power. 

To avoid numerical simulation, different simplified models of light trans¬ 
mission and reflection through ID gratings under sub-wavelength (w < 0.5A /n g , 
where n g is the dielectric constant of the slits) condition have been developed 
(l, 3, 5|. In developing these models, the fundamental waveguide mode is the 
only propagating mode inside the slit/cavity is assumed. In addition, these 
models also assume that the grating ridge and/or the slit walls are made of 
Perfect Electrical Conductor (PEC) [1, 3, 5]. For most of the metals the 
PEC assumption holds true in the deep infrared and microwave regions of 
the electromagnetic spectrum and that is where the previous models produce 


reliable results. In the visi 
become significantly lossy 


e band of the electromagnetic spectrum, metals 


, |17] and the PEC assumption along with the ex¬ 
isting transmission and reflection models of light fail. In this context, a model 
which does not assume the PEC assumption is necessary to analyze ID grat¬ 
ing structure in the visible band of the electromagnetic spectrum. Further, 
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a model of ID metallo-dielectric gratings which can provide explanation of 
light transmission and reflection in terms of grating geometrical parameters 
is also relevant considering the limitations of the existing models. In this 
article, we aim to develop a model of ID metallo-dielectric gratings of real 
metals which complements the shortcomings of the current semi-analytical 
(involves infinite sums) models. In developing this model, w < A/(2 n g ) is 
assumed. 



Figure 1: ID Lamellar Grating 


2. One mode model of ID transmission grating 

Let us assume a TM polarized electromagnetic wave H incy = exp{—icut}* 
exp {mifc 0 (sin 6tr — cos#,z)} of wavelength A is incident upon the grating 
structure shown in Fig. |T] at an incident angle 9, where n\ is the refrac¬ 
tive index of the incident medium (z > 0) and k Q = 2tt/X is the free space 
wave number. n g and n r are the refractive indexes of the grooves/slits and 
the grating ridges respectively. The dielectric constants of the slits and ridges 
can be written as e g = n 2 and e r = n 2 correspondingly. A part of the in¬ 
cident energy is reflected back to the incident medium as diffraction orders 
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whereas another part emerges as transmitted diffraction orders in the sub¬ 
strate (—h < z ) of refractive index n s . Using Rayleigh expansion the reflected 
and transmitted diffraction orders can be expressed as - 


jjr tt 

y 11 inc,y 


+ Rm exp [ik 0 {p m x + {n\ - 7m ) 1,2z }] 

m 

H l = Tm exp i ik oHnX - (n 2 s - 7m) 1/2 -}] 


( 1 ) 


where 7 m = ri\ sind + mX/P and m is an integer. R rn and T m are the 
reflected and transmitted magnetic field amplitudes fl^ |. 

The transmitted energy passes through the slits from the incident medium 
to the substrate as waveguide modes. Based upon A, w, and P, different 
waveguide modes are excited inside the slits. Mode profiles and the x— 
and components of wave vectors of different waveguide modes can be 
found by solving Maxwell’s equations in the grating region {—h < z < 0) 


in combination with the stratified medium theory 


14], The x— and z— 


components of various waveguide modes can be found from Eq. 


cos(fco^i-Psind) — cos (/3tP) cos (afP) H—sin (/3tP) sin (afP) = 0 (2) 

2 p e r a 

where t — (P — w)/P, f — 1 — t — w/P, and a = k 0 e g — n 2 E ^ and (3 = 
ko\j e r ~ n2 Eff are x ~ components of a waveguide mode in the slit/groove 
and ridge materials respectively. The z— component of a waveguide mode 
is defined as k z = riEffko. For metallic ridges, the dielectric constant is 
given by e r = n 2 = (rj + in) 2 , where 7 and k are the real and imaginary 
components of the refractive index, respectively. As can be understood from 
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Eq. ((2D that there are infinite number of solutions of this equation and each 
of these solutions is called a waveguide mode, ffowever, when w < A/2 with 
n g = 1, only the fundamental waveguide mode is propagating inside a slit 


cavity 


1-3]. In this case, given that \nj\ >> riEff of the fundamental 


mode can be expressed as Eq. 


18) 


n 


Eff 


= 1 + i- 


\ 


7 Twn r 


(3) 


Considering Eq. (J3]) , a and {3 can be expressed as k 0 ^/—i\/(nwn r ) and 
k 0 n r respectively. Mode profile of the fundamental waveguide mode can be 


found from Eq. 


14]- 


v , \ tt r ( tP \i , - v ^ e 9 k o ■ r i tP \} 

X(x) = U\ cos [a(x -—)J + j --—sm [a(x -—)J, 

la l 

\ of , , - v o£rk 0 . r . tP 

= cos[/3(x+—)]+j—sm[^(a; +—)], 

where Vo, Vi, Ui, M and N are give by 

V 0 = [exp(jk 0 P sin 6) — M\/N (5) 

Vi = j-r 1 — sin ((3tP) + Vo cos ((3tP) (6) 

k 0 e r 

U x = cos (ptP) + j k ° V ° 6r sin (ptP) (7) 

M = cos (/3tP) cos ( afP )-— sin (/3tP) sin ( afP ) (8) 

6 r OL 

N = jA;o[cos (PtP) sin (afP)/a + e r sin (ptP) cos (afP)/P] (9) 

Given that | a/— 247ru;/(n r A)| << 1 , | cos yJ—iAnw/(n r A)| m 1 , and 

| sin \J—iAnw/(n r A)| ~ 0, M and N can be expressed as below -. 


Y S M < t 1 - 

tP . , 

1^1 < y ( 4 ) 
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M 

N 


cos {k 0 (P — w)n r } cos (\l —i ——) 

u A n r 


in r sin {k 0 (P — w)n r \ cos (\ — i ——) 

A n r 


. Anw . 


( 10 ) 

(ii) 


Further, provided that | exp {ifco-Psin #}| = 1 and \N\ » 1, Vo, V\ and 
Ui can be expressed as below 


V, 

U { 

V, 


0 

.cot {k 0 (P — w)n r } 


n r 


( 12 ) 

(13) 

(14) 


Substituting M, N, Vo, Vi, U\, a and (3 in Eq. (J3J, X(x) corresponding 
to the fundamental mode can be expressed as below 


X(x) ~ cos{fc 0 ((r + 


X{x) 


w P 

— ) n r} ~ cot {k 0 (P — w)n r } sin {k 0 (x H- 


tP tP 

- < x < — 

2 ~ ~ 2 

tP . . t. 

T- W -(1 ~2 )p 


W , 


»,) 


(15) 


Provided that the fundamental waveguide mode is the only surviving 
mode, the magnetic field inside the grating region (— h < z < 0) can be 
expressed as 


H gy = X(x)[Aexp (in E ffk 0 z) + B exp (—in E ffk 0 z)] (16) 

where A and B are the co-efficients of the counter propagating fundamen¬ 
tal mode inside the slit cavity. After applying boundary conditions at z = 0 
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and z = —h interfaces for the magnetic and electric fields in combination 


with Maxwell’s equations, one gets relations involving A, B, R m and T m as 
follows 


{A + B)X m — R rn = 5 m o (17) 

[A — B)Q m — Tl m R m = — fI m 5 m o (18) 

{AT, + BT 1 )X m — A m T m = 0 (19) 

{AT — BT 1 )Q m + A m £ m T m = 0 (20) 


where T = exp{—ihk z ) and T 1 = exp{ihk z ) are constants. k z = riEffko 
is the ^-component of the wave vector corresponding to the fundamental 
waveguide mode. fl m = k Q y/(n\ - 7 ^), £ m = k 0 ^{n 2 s - 7 ^)/^, and A m = 
exp {ikohsjn 2 s - 7 ^). X m = (1/P) X{x) exp {-ik 0 'y m x)dx, Vt m = 

{k z /P) f^ P ^ P f^(ic) exp {—ik 0 ^/ m x)dx, and h2(a;) = X{x)/e. For — tP < x < 
tP , e is equal to e r whereas for tP < x < (1 — |)P, e is n 2 g = 1. 

After some simple algebraic manipulation of equations (TTSj)- ([20]) and as¬ 
suming ri\ — 1, T m and R m can be written as follows. 


T = 

1 m 


2 Smo exp (- ik 0 hyjn 2 s - 7 ^) 


i 1 + i cos o'*) - + AT-hJ sin » ft) 


( 21 ) 


Rrn 


2 h m0 [cos ( n E ffk Q h ) - sin (n E// fc 0 h)] 

t 1 + 77 ] C ° S ( nE B k ^ ~ + £r 7i-7,iJ Shl 

—^mO (22) 


By expressing ng// as ri E ^ -\-in E ff, and given that sinh (n^Poh) << 1 






















(since n' E jj « 1) and | —|^| << 1, for m = 0, T 0 and i? 0 can be simplified 
to - 


T n = 


2 exp (— ik 0 hn 2 s — sin 2 6) 


[(! + h VlS?) cos ( n EfM) - i €rV ,^f 9 sin (n' Eff koh)} 


x 


cosh (n E fj-koh) 


(23) 


Rn — 


2[ cos ( n Eff k 0 h ) - 6 sin (n' Eff k 0 h)\ 


( X + 


n " sm2g ) cos ( n' Eff k 0 h ) - ° sin ( n Eff k 0 h ) 


1 (24) 


l-sin^e J^VEff 

It can be seen from Eqs. CD and 
analytical models 


n E ff \--Eff 

that unlike the existing semi- 


L|, [3, 5] of light transmission and reflection through the 
ID grating structure, the models presented here are analytical and hence 
the physics of light transmission and reflection can be easily understood. In 
the next section different examples are considered and other benefits of the 
proposed models are discussed. 


3. Results and Discussion 


To justify the validity of the transmission and reflection models developed 


in the last section, geometrical grating parameters from Ref. 


19] are consid¬ 


ered. Specifically, P = 75 nm, w — 20 nm, h = 198 nm, n s — 1, 9 — 0° and 
the grating is made of gold are assumed. These parameters are supposed to 
produce a Fabry-Perot resonance peak at 750 nm and no transmission peak 
due to the Wood-Rayleigh Anomaly (WRA) or the surface plasmon resonance 
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Figure 2: Zeroth order (a) transmission, and (b) transmission and reflection spectra of a 
gold grating for P = 75 nm, w = 20 nm, h = 198 nm, n s = 1, and 9 = 0°. 
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is expected in the visible band. Fig. [2(a) shows the zeroth order transmis¬ 
sion spectrum corresponding to these grating geometrical parameters. In this 
hgure, in addition to the data available from Eq. (l23]h results obtained from 
the exact modal analysis Q with thirteen modes and the analyt cal model 
developed by Garcia-Vidal et ah (5 have also been included. One can see 
that, except the amplitude mismatch, the analytical model developed here 
(Eq (1231) ) fits very well with the EMA results. In contrast, as mentioned 
earlier, the semi-analytical model developed by Garcia-Vidal and co-workers 
5] completely fails to predict the resonance. However, if one closely looks at 
the denominator of Eq (l23lh then it is evident that when the FB resonance 
condition is met ( n E ^k 0 h = jn with j = 1,2,3,.. 2]), cos (n E ^k 0 h) = 1 
and sin (n E ^k 0 h) = 0. The consequence of the FP resonance on the denom¬ 
inator is that it becomes smaller compared to non-resonance scenarios. This 
ensures that To is maximized and a transmission peak appears. Fig. [2](b) 
confirms that at the transmission peak, reflection is reduced as expected. 
Fig. [3] shows the dependence of the zeroth order transmission on the angle 
of incidence. It is known that when P « A (in our case P = 75 nm and 


A = 750 nm), the dependence of the FP resonance on 6 is minimal 


[lit 


2Q|. 


This is confirmed from the results of Fig. [3J It can also be understood from 
Eq. (l23]h Elaborately, irrespective of the incidence angle, the absolute value 
of the numerator of Eq. (123|) is always 2. Additionally, when the FP reso¬ 
nance condition is satisfied and n s = 1, the 1 st and the 2 nd terms within the 
square bracket of the denominator of Eq. (1231) become independent of 6 and 
zero, respectively. Upon these circumstances, Eq. (123|) simplifies to Eq. (l25]h 
From this last equation it can be seen that the incidence angle has no impact 
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A, (nm) 


Figure 3: Zeroth order transmission at different incidence angle for the same grating 
parameters of Fig. [2j 


on the FP resonance. Further, front Eq. (1231) it can be understood that a 
peak related to the FP resonance in the extraordinary optical transmission 
through a ID metallo-dielectric grating is independent of the grating period 
(P). This has been experimentally verified by Pang and coworkers 2l|. Also, 
it is evident from Eq. (|25|) that as the grating thickness (h) increases (for 
further discussion on the impacts of h on transmission, see below), the in¬ 
tensity of a peak related to the FP resonance decreases. At this point it is 
relevant to mention that these are the kind of insights that our model can 
provide while other existing models of ID metallo-dielectric gratings fail to 
deliver. 


To 


1 

cosh (n'zffk 0 h) 


(25) 


The role of the slit width on the transmission intensity can be appreciated 
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through the influence of n" E ff on Eq. (125]) . It is known that varies 
inversely with the slit width {w , see Eq. (J3D) or n' E ff increases as w decreases 
and vice versa. The result of the inverse relationship between w and tl E jj is 
that as the slit width shrinks cosh ( n^^k^h ) in Eq. (123]) increases and hence 
the transmission intensity diminishes as intuitively expected. 

To understand the effects of the grating thickness (h) on the overall trans¬ 
mission intensity related to the Fabry-Perot resonance, let us go back to Eq. 
(121]) with 6 = 0°, n s = 1 and m — 0. Neglecting terms with —in the 
denominator, one gets - 


T n = 


2 exp(— ik 0 h ) 


2 cos (■ n Eff k 0 h ) - sin (■ n Eff k 0 h ) 

T 0 = 2 exp (—ik Q h)/ 


[(2 cos (n E fjk 0 h) — i 


n Eff 

6 r 


■ sin 


n Eff 
£ r 


(n Eff k 0 h )) cosh ( n" Eff k 0 h ) - 


(2 sin (n E ffkoh) +i -cos (n E ffkoh)) sinh (n E ffkoh)] (26) 

It is clear that as h increases both cosh {n' E ^jkah) and sinh (n E ^koh) in 
the denominator of Eq. (1261) rise rapidly and as a consequence T 0 is expected 
to reduce significantly. In other words, if one uses a thicker grating to attain 
a FP resonance at a desired A then the associated transmission peak will 
be weaker compared to the grating which has a smaller h. This has been 
experimentally confirmed by other research groups P . I22I ]. Lastly, to un¬ 
derstand the role of the grating material (specifically the ridge material) on 
the transmission through a ID grating, inspecting the denominator of Eq. 
(126]) can again be useful. Elaborately, for a fixed w and P, as e r increases, 
the denominator of Eq. (126]) increases and consequently T 0 is supposed to 


13 











Figure 4: Zeroth order transmission for P = 418 nm, w = 50 nm, h = 209 nm, n s = 1.52, 
0 = 0° and the ridge metal is gold. 


decrease proportionately. Alternatively, it can be concluded that if one con¬ 
siders two metals of different dielectric constants to get a FP resonance with 
the minimum possible grating thicknesses at a desired A (a metal with a 
higher n r requires a thicker grating to fulfil the FP condition, see equations 
(]3l) ) then the metal which has a higher |e r | (in particular a metal which has 
the higher imaginary dielectric constant) will produce a weaker transmission 
peak compared to the metal which has a smaller \n r \. Again this has been 


verified 


9], 


Before concluding, one more example is considered. In this example, we 
consider P = 418 nm , w = 50 nm, h = 209 nm, and a gold grating 23] 


on a glass substrate (n s = 1.52). These grating geometrical parameters are 
supposed to produce a transmission peak at around 635 nm. Equivalently, 
the grating should work as a color filter by transmitting the red light while 
blocking all other colors 24], Corresponding to these parameters, a Wood- 
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Rayleigh anomaly and a FP resonance should occur at 635 nm. Fig. H] shows 
the zeroth order transmission spectra obtained using the model ( Eq. ([23]) ) 
developed in this article, the model of Lalanne et al. jlij] and the Rigorous 
Coupled Wave Analysis (RCWA). It can be observed that, like the RCWA, 
the simple analytical model developed in this article, despite some ampli¬ 
tude mismatch, can predict the wavelength and the profile of the resonance 


f a 


while the model of Ref. [3!l completely fails. According to us the main rea¬ 
son behind the failure of Lalanne and co-workers’ model of ID grating is 
the PEC assumption that they have made in deriving their model. As men¬ 
tioned earlier, the PEC assumption is valid deep inside the infrared and the 
microwave regions of the electromagnetic spectrum and that is where this 
model matches very well with the RCWA results (such as those of Ref. 6]). 
In the visible band of the electromagnetic spectrum, the perfect electrical 
conductor assumption is invalid for metals like gold, silver and copper. Con¬ 
sequently, results produced using the model of Lalanne et al. are inaccurate. 
Lastly, it is noticeable from Fig. [4] that the line-width of the peak of the ana¬ 
lytical model developed here is broader than that obtained from the RCWA. 
Such a broader peak is normally a characteristic of the FP resonance. Since 
the model developed here is insensitive to the Wood-Raylcigh anomaly, the 
effects of the WRA is not reflected in the spectrum obtained using Eq. (l23j) . 
Consequently, a mismatch between the result generated using the RCWA, 
that includes the influences of both WRA and FP mechanisms, and the data 
available from Eq. (l23j) is expected. 


15 




4. Summary 


In conclusions, an analytical model of light transmission and reflection 
through ID metallo-dielectric gratings based upon the exact modal analysis 
has been developed. This model can accurately predict a transmission peak 
due to the Fabry-Perot resonance. It has also been theoretically proved 
that transmission intensity of a peak related to the Fabry-Perot resonance 
decreases with the increasing grating thickness and the refractive index of the 
ridge material. Further, when the grating period is much smaller compared to 
the incident wavelength, the model presented in this article can capture the 
variation of transmission intensities as the incident angle changes. Limitation 
of this model is that it can not stipulate a transmission peak due to the 
Wood-Rayleigh anomaly. 
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